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SZ effect or Not? - Detecting most galaxy clusters' main 
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ABSTRACT 

Galaxy clusters are the most massive objects in the universe, and they comprise a high 
temperature intraclustcr medium of about 10 T K, believed to offer a main foreground 
effect for the CMB data with thermal Sunyaev-Zel'dovich (SZ) effect. This assumption 
has been confirmed with SZ signal detection in hundreds of clusters, but comparing 
with the huge numbers of clusters within optical selected samples from SDSS data, 
this only accounts for a few percent. Here we introduce a model-independent new 
method to confirm the assumption that galaxy clusters offer the thermal SZ signal as 
their main foreground effect. For the WMAP 7year data, we classified data pixels as 
"to be" or "not to be" affected by the sample clusters, with a parameter of its nearest 
neighbor cluster's angular distance. By comparing the statistical results of these two 
kinds of pixels, we can see how the sample clusters affect the CMB data directly. We 
find that Planck-ESZ sample and the Xray samplcs(^ f 2 clusters) can lead to obvious 
temperature depression in WMAP 7year data, this confirms the SZ effect prediction. 
However, each optical selected sample (> 10 4 clusters), shows an opposite result: the 
mean temperature rises to about 10 uK. The unexpected qualitative scenario implies 
that the main foreground effect of most clusters is NOT always the expected SZ 
effect. This is maybe the reason why the SZ signal detection result is lower than what 
is expected by the model. 
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1 INTRODUCTION 

As the most massive self-gravitating systems in the cos- 
mos, galaxy clusters can make significantly contribute to the 
mensuration of precision cosmology and during their forma- 
tion and evolution various effects can be analyzed statisti- 
cally. One major effect is t he Sunyaev-Zel'dovich (SZ) effect 
ISunvaev-Zel'dovich 19721 ) : the CMB photons can have in- 
verse Compton scattering with high energy electrons in in- 
tracluster medium when passing through clusters. The ther- 
mal SZ effect is considered to be the most remarkable effect, 
as it can increase photons' energy statistically and notice- 
ably distort the CMB spectrum. The thermal SZ effect de- 
creases intensity in low frequencies ( like Q, V, W band of 
WMAP) and increases intensity in high frequencies. 

This predictable distortion of CMB data is expected as 
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the marked signal of clusters. After the SZ signal of res- 
idential clusters was confirmed, the blind sky surveys us- 
ing SZ e ffect have continue d over t he last decade, includ - 
ing SPT (I Arthur et al. 20061 ). ACT (|Carlstrom et al. 201 j ). 
PLANCK~( |pianck 20061 ) and so on. A number of high red- 
shift clusters are expected to be found via the blind SZ sky 
surveys using feat ures of SZ signal and it s insensitive prop- 
erty on red shift (|Carlstrom. et al. 20021 ). This provides a 
solid foundation for the mensuration of precision cosmology. 

Here we notice a basic assumption is behind the ex- 
pectations of these ongoing projects: the main foreground 
effect of most galaxy clusters to the CMB data, should be 
the expected thermal SZ effect. Though people have already 
detected and confirmed hundreds of clusters' SZ signal, we 
cannot assume that this is the case for all clusters. Compar- 
ing with the huge cluster numbers in optical selected sam- 
ples from SDSS data, the basic assumption has only been 
validated in a few percent o f clusters (and it is hard to pr o- 
vide their selection function jPlanck Collaboration. 2011al )). 
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Considering its unknowns can have unpredictable or unfore- 
seen impacts on understanding or applying the results of 
SZ effect galaxy cluster survey, therefore we need a direct 
detection of most clusters to confirm this assumption. 

At the same time, the cosmological analy- 
sis rnethojfs_sucli__as__^hister number count" et 
al. (|Gilbert Holder et al. 200ll ; iManera fc Mota 20061 ') 

expect a complete catalogue of setting conditions. It's 
important to make sure that the SZ effect signal blind 
survey of galaxy clusters will not miss cluster samples. 
However, some studies have already found the observed SZ 
effect signal to be (es pecially for optical selected clusters) 
not strong enough (|Bielbv et al. 20071; iLieu et al. 2009; 
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Planck Collaboration 
2012I ; 



2011al; 



Diego et al. 200S 



Draper et al. 



Neelim a Schaa l et al. 2012f)fthough 
debate ( Planck Collaboration. 2011bl : 
Melin 20 111 )') .which rm- 



still in 

lAfshordi. Niavesh et al. 200?! : 
plies the existence of other considerable foreground effects 
of clusters will contaminate the SZ effect signal. The 
existence and negligibility of SZ effect signals is becoming a 
noteworthy debate. We notice that some traditional analysis 
methods are model-dependent and the free parameters 
can lead to uncertainty in the debate. It is necessary 
to introduce a model-independent new method that will 
ensure more reliable conclusions can be drawn. 



2 METHODS 

For the main sample galaxy cluster we used typical size 
of five Mpc/h at z ~ 0.5 distance, its angular size (for 
ACDM model) is about six arcmin. This is round one pixel 
size of a Healpix data (iV S jd e =512). So these objects ap- 
pear as "points source" and will quickly "disappear" in the 
background noise when convolved with a beam as wide as 
WMAP's beams. It is hard to detect and confirm each single 
cluster's SZ signal. However, their statistical effect will not 
disappear. One can take help from the viewpoint of CMB 
data pixels, just take each pixel as a probe. To study the 
foreground effect of one galaxy cluster in an ideal isotropy 
CMB, a simple method is used to compare the probe data 
(temperature data of this pixel) of angular regions affected 
by the cluster or not. For real CMB data, the fluctuation 
temperature of each pixel can be taken as the detector's 
another Gauss distribution error. Considering the different 
properties of noise signals and the SZ signal, one can use 
statistical method to compare the mean probe data of an- 
gular regions of "to be" or "not to be" regions affected by 
the sample clusters. The noise signal will have similar effects 
on these two kinds of pixels but the thermal SZ signal will 
only depress the temperature of "to be" affected pixels. 

The preconditions here are two points: First, we are 
able to assort these two angular regions; Second, each region 
includes enough data pixels to minimize the statistical error. 

We find a continuous parameter di is competent for 
such taxology. For each CMB data pixel, di is defined as 
the angular distance of this pixel to its nearest neighbor 
galaxy cluster (of the cluster sample used). Comparing the 
traditional "stacking method" can help us understand the 
parameter di. The stacking method selects each cluster as 
the origin and bin pixels by circular rings around it (with 
different angular distances), and stacking annulus of all clus- 
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Figure 1. Pixel region round a group of cluster within a bin value 

of 9' < di < 13'. Each trigon correspond to the position of one 
cluster in GMBCG sample, and each gray disc means the central 
position of one pixel in WMAP data. For sparse cluster samplefthe 
up-right one), we can see our method will simply retrogress to the 
stacking method. 



ters to do analysis. Considering that galaxy clusters tend 
to swarm together, the temperature signal of pixels within 
annulus round one clusters can be seriously affected by its 
neighbor clusters. In Fig[TJ we show the pixel region around 
a group of clusters within a bin value of di. Our method can 
be comprehended as changing the circle annulus around each 
cluster with curved loops around the collection of sample 
clusters, and redefining the angular parameter as d\. This 
means the angular distance to the whole cluster sample is 
used rather than with each cluster. 

As clusters mainly affect their local angular region, pix- 
els with di parameter large enough can properly represent 
the "not to be" affected angular regions unlike the stacking 
method. One other technical merit here is the impartial use 
of pixel data, as the stacking method applies pixels affected 
by neighbor clusters more times. When the cluster sample is 
sparse in the sky, our method simply reverts to the stacking 
method. 

H ere we illustrate the result s with optical selected GM- 
BCG (jjiangang Hao et al. 20101 ) galaxy cluster sample. In 
Fig(2] we show the distribution function of pixels' di pa- 
rameter. The figure shows the effective statistical region is 
within < di < 40'. By comparing the possible angular size 
that one cluster might affect pixel data (mainly the beam 
size), one can see both the cluster affected region and un- 
affected region that contain enough pixels for effective sta- 
tistical analysis. In contrast, due to the angular resolution 
of WMAP data, this is difficult to do with SDSS galaxy 
samples. 

We can do analysis with the < T >-< d\ > curve (here- 
after TD curve), by taking each pixel as a probe, comparing 
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Figure 2. The distribution function of pixels' d± parameter. The 

pixels correspond to the CMB data within the main angular region 
of GMBCG sample. We select the 50580 clusters of the sample in 
the main survey area of SDSS DR-7 and take 503740 WMAP data 
pixels within the area as our statistical CMB angular region. In order 
to avoid some unwanted edge effect, we have already dropped pixels 
close to the edge within 70'. The blue lines show the beam size of 
WMAP data in Q, V, W band, respectively. And the percentages are 
the proportions of pixels whose d\ value are within these half beam 
size. 



the mean temperature < T > of pixels binned with differ- 
ent di value. The merits of this method will be discussed 
in detail in another paper. Here we emphasize two points in 
physics: 

First, the two sides of the TD curve represent different 
cases of pixels being affected or not (by the cluster sample 
used). The main foreground effects of galaxy clusters that 
we are interested in (such as SZ effect and radio emission) 
have the property of "angular localization", which means 
they only affect the angular region they appear in (within 
several arc-minutes for most clusters). Since the beam angle 
of WMAP data in Q,V,W band and ILC data ranges from 
13' to about 30', it is safe to say pixels of d\ — > 40' represent 
angular regions "surely unaffected" and the di ~ 0' pixels 
represent " affected regions" by the cluster samples. By com- 
paring the mean temperature of these two cases, we can see 
how galaxy clusters affect CMB data directly. 

Second, some background or foreground unrelated ob- 
jects outside the cluster sample might also affect the CMB 
data, but statistically they'll change the two side of the TD 
curve in the same way. This can be tested by simulation. So, 
if a reliable difference on the two side is confirmed, it should 
be the effect caused by the cluster sample itself. 

Now we can see the qualitative scenario about how the 
cluster sample affect the TD curve: When thermal SZ effect 
being the main foreground effect, the mean temperature of 
WMAP data in Q, V, W bands should noticeably decline for 
cluster regions(di — > 0'), and the large di side should remain 
unaffected. Then, the TD curve should rise from d\ = 0' 
to di — 40'. In contrast, if radio emission was the main 
foreground effect, the low di side would be driven up and 
the TD curve would be decline. 
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Figure 3. Comparing the mean temperature of pixels with differ- 
ent di parameter. The figure show typical TD curves of different 
CMB data - galaxy cluster samples. The temperature distribution 
function of pixels within each di bin (3') can be well described 
with a Gauss function, and we calculate the error bar of < T > 

/y^(T ! -<T>) 2 

with a = y ' n(N—i) ' s tne P' xe ' number within the d\ 

bin). The curves can be described with an empirical line function: 

< T >= A* < di > +Ti (fitting with weight \/~N). The tradi- 
tional stack method is somehow our special case when the cluster 
sample is sparsity and the neighbor-cluster effect can be neglected, 
this is the case of Planck-189 sample. For the Planck-189 sample, 
we select WMAP data pixels round each cluster within 55' as our 
statistical region (bin-width is 5.5'). Since only the slope A is really 
interested , the WMAP(ilc)-Planck(189) cluster TD curve is trans- 
lated as < T > — 13mA' here, to show more clear contrast with same 

< T > in unaffected region. 



3 RESULTS 

In Fig[3] we illu strate the TD curve of optical selected 
GMBCG sample ijjiangang Hao et al. 20101) and Planck SZ 
effect selected 189 clusters ( Planck Collaboration. 2011ai ) 
(Planck-189). For Planck-189 sample, the TD curve line is 
obviously rising, which means < T > is much lower when 
d\ — > 0' than that in large di region, thus well confirming 
the SZ effect prediction. However, when the foreground is 
changed to the 50,580 clusters of GMBCG sample, we get 
an unexpected opposite result: the TD curve becomes a vis- 
ible downward curve, which means the mean temperature 
behind clusters has got an increment up to about lOuK. 
This result is similar for WMAP data in Q,V,W bands and 
also the ilc data, with small error margins. (While, the TD 
curve becomes a nearly level line when we use the CM B data 
eliminated by another team |Delabrouille et al. 20091 ) using 
needlet method.) 

As a post hoc examination, we also calculated the TD 
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Table 1. Fitting results of different samples. 



curves using the simulated CMB data of WMAP, and also 
the simulated random distributing cluster sample. For these, 
no causal relationship data were found, Fig(4]shows their TD 
curves are common level curves as expected. 

Fig'Oalso shows us that the TD curves of different sam- 
ples are typically approximately straight lines, so we can 
describe them with an empirical line function: 

< T >= A* < di > +Ti 

If we are not using one all-sky cluster sample, the param- 
eter Ti can be influenced by the cluster sample region and 
also the CMB large scale fluctuation. The valuable param- 
eter is the A slope. For the samples we used, the value 
ATa = —^4*40' can roughly show us the difference of < T > 
between affected and unaffected regions, which underlies 
how these cluster samples affect the CMB data. ATa < 
correspond the result of thermal SZ effect and ATa > 
relate to opposite effect like radio emissions. 

In tabl^T] we show the value of ATa when setting 
different cl uster samples as foregro und of WMAP 7- 
year data (|Gilbert Holder et al. 201 ll ) (W band and ilc 
data) . Here we can see the ATa values fall explicitly 
in two situations: for the SZ effec t selected and Xray se - 
lected cluster samples (ACT jMarriage. Tobias 20lj ) , 

Planck-189 (iPlanck Collaboration. 2011al),xcs 3 

(|Mehrtens. Nicola et al. 201ll ).mcxc (|Piffaretti et al. 201ll )) 
, ATa is obviously negative, confirming the SZ effect 
image; yet for each optical selec ted clus ter samples (GM - 
BCG Jjiangang Hao et al. 20l"oT ). Wen (|Wen et al. 20121) , 
maxBCG~i |Koester et al. 20071 ) ), the ATa value is signifi- 
cantly positive. 



4 DISCUSSION AND CONCLUSIONS 

Before quantitative study about ATa value in more detail, 
we focus on the explicit ATa > property of each optical 
selected cluster sample. Its value is about zero when we use 
a random distributed foreground sample in Fig(3] It is inter- 
esting to understand w hy ATa changes in Fig|3] The ATa < 
case of Planck-189 (|Planck Collaboration. 201 lal ) corre- 
sponds to thermal SZ signals directly, then the unexpected 
ATa > case of GMBCG sample with their main fore- 
ground effect on WMAP data is NOT the expected SZ effect. 
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Figure 4. The TD curves of testing data of CMB and clusters. 

This figure includs: 1, randomized cluster sample (of GMBCG); 2, 
simulated CMB data in Q,V,W band; 3, reversal of CMB ilc data 
rl80(ra -> ra+ 180°), n2s(dec -s> -dec). For all the CMB data, we 
select the 503740 pixels within GMBCG's main region. (Since it's not 
using all-sky pixels, we are common to see their mean temperature 
is not zero.) 



Such a result implies the existence of serious radio emission 
like effect in these optical select clusters, and this scenario 
clarifies w hy the SZ effect signal is appa r ently weaker than 
expected JPlanck Collaboration. 2011a : Diego et al. 20031 : 



iDraper et al. 20121 : iNeelima Sehgal et al. 20121 ) 



• In angular regions of galaxy clusters, there exists an 
opposite contamination foreground effect something like ra- 
dio emission that can affect the SZ effect signal detection. 
Such an opposite signal in most clusters is high enough that 
it can cover the SZ effect signal of the cluster, and being 
the main foreground effect. It should not be neglected when 
doing CMB signal analysis(soever they are surely clusters or 
not). 

• In distance, this contamination should come from the 
cluster itself. Because if it affects on the line of sight before 
or after the clusters, such as the star burst galaxy at high 
redshift background, statistically there should also be the 
same effect at the non-cluster regions, and not result in a 
temperature change. 

• In spectrum, such "emission components" show simi- 
lar effects on Q,V,W band, a little higher at high frequency. 

In conclusion, our model-independent method shows 
the main foreground effect of most galaxy clusters directly. 
The results of known SZ signal selected clusters and Xray 
observed clusters confirm the traditional thermal SZ result. 
Unexpectedly, however, most clusters in optical selected 
samples' thermal SZ signal are contaminated (even covered) 
by something like radio emission. This may be the reason 
why the SZ signal detection result is lower than model expec- 
tations. This may also help us understand the cooling flow 
problem of galaxy clusters: if the unknown effect is surely 
radio emission, then it can imply the mean radio intensity 
in a special way, for these clusters are able to offer radio 
emissions (on average) at least covering their thermal SZ 
signal. 
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